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ABSTRACT 

We use the Wilkinson Microwave Anisotropy Probe 7-year data (WMAP7) to further probe point source detection 
technique in the sky maps of the cosmic microwave background (CMB) radiation. The method by Tegmark et al. for 
foreground reduced maps and the Kolmogorov parameter as the descriptor are adopted for the analysis of WMAP 
satellite CMB temperature data. Part of the detected points coincide with point sources already revealed by other 
methods. However, we have also found 2 source candidates for which still no counterparts are known, and identified 7 
point sources listed in Planck Early Release Compact Source Catalogue as high reliability sources. 
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1. Introduction 



The cosmic microwave background (CMB) maps have been 
examined for detecting point sources as potential fore- 
grounds, which could be cosmological or Milky Way objects 
emitting thermally or non-thermally. A catalog is made 
available on the NASA WMAP team web page. Various 
methods, including the wavelets and needlets, have been 
used for t he detection of point sources in pixel i zed sk y maps 

(|2011|)). Most 



(see, e.g., Scodeller et al. (2012 
of the detected sources in CM 
radio sources, quasars, blazars 



Batista et al. 

maps coincide with known 
although some sources a re 
still unidentified (Wright et al.||2009[ parosik et al.||2011 ). 



Kolmogorov stochasticity pa rameter (KSP) has been 



already involved for such aims (Gurzadyan et al. 2010) 



some of the sources revealed by that method initially had 
no counterparts, however their counterparts have been de- 
tected later by the Fermi satellite as gamma sources and 
were included in its first source catalog 1FGSC. In the 
present study we continue the application of KSP technique 
to detect point sources in WMAP7 data set, and for that 
aim we borrow the maps cleaned by ( jTegmark et al.| [2003) 
method. While applying Kolmogorov method we analyze 
the mask issue, and we study the cross-correlations between 
the power spectra of the these two types of maps, since the 
K-maps can carry information about the matter distribu- 
tion in the Universe, particularly on voids ( Gurzadyan et al. 
2009bl. 



This paper is organized as follows. First we introduce 
briefly Kolmogorov method, then we construct cleaned K- 
maps using Tegmark et al. method, modified for such maps. 
We use the resulting maps to reveal the point sources, and 
to discuss the cross power spectrum between the CMB tem- 
perature and Kolmogorov maps (K-maps). 



2. Kolmogorov method and CMB maps 



In his original work, Kolmogorov ( 1933 ) introduced a 
method for determining whether the given random num- 
ber sequence X{,i=l,..,N, ordered in an increasing 
way, obe ys to a given statistic or not (Arnold 2008a|b|c 



2009a b). For such a purpose the so-called empirical cumu 



iative distribution function F^(x) is calculated, 



F N (x) 




x < Xi 
Xi < X, k 
Xn < x 



1,2,. 
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where k is the number of elements which obey to the rela- 
tion Xi < x. Having assumed a particular theoretical cu- 
mulative distribution function (CDF) F(x), the parameter 
Xn is easily calculated, 



A 



A? 



In sup \Fff(x) — F(x) 



Kolmogorov proved that in the limit N — > oo, Xn, which is 
a random variable, has a cumulative distribution function 
$(A) reading as 



$(A) 
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The function $(A) is actually also a particular value of 



thet a functions, since $(A) = ^4( 0, e 2A ) = $ 3 (0, 
(see |Abramowitz fc Stegun||1970[ 16.27). 



-2A' 
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3. Cleaned CMB K-map 

3.1. WMAP7 W CMB and K-maps 

It is well known that the CMB sky has an nearly Gaussian 
distribution. So, the Gaussian distribution is used in 
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Kolmogorov method to construct a K-map. For every com- 
pact region of the CMB map containing 256 neighboring 
temperature pixels, one corresponding pixel value of the 
K-map is obtained. For WMAP7 W band C MB map with 

^side 

( |2005 )), we obtain a K-map with n.: 1A i, 



512 resolution parameter (details in Gorski et al. 

32 parameter. 



rhis means that CMB map has n map = 12rig idc = 3145728 
pixels and the K-map has n map = 12288 ones. This is due 
to the fact that KSP is a statistical parameter. Then the 
KSP distribution maps over the whole sky can be obtained 
for WMAP7 data. It is seen that the Galactic disk region 
has higher and saturated KSP values, which indicates that 
it has a non-Gaussian distribution, distinct from the CMB. 
Also a lot of pixels have a high value of KSP, but most of 
them, as it will be shown below, arc due to instrumental 
and other types of noise, also of non-Gaussian nature. 

Throughout the analysis we mainly use the HEALPIX 
program package for calculation of the spherical harmonic 
coefficients and to reconstruct the maps |G6rski et al.| 



(2005). For some tedious manipulating procedures with 
the spherical harmonic coefficients ai m (see Eq. (J6| below) 
we us e the GLESP program package |Doros hkcvich et al.| 
(2005). For maps with low resolution (in our case for K- 
map, riside = 32,n map = 12288) the ai m calculation and 
map reconstruction the use of HEALPIX program package 
led to the calculation accuracy problems. To test whcthcr 
the calculation error is small or not for a low resolution 
map, we construct a unity HEALPIX map with n si( i c = 32, 
and another one with n s id c = 512, and run them through 
the ai m calculation and map reconstruction procedures. 
The result, multiplied by 1000, is given in Fig. [2] For both 
cases this procedure adds some non-isotropic noise which 
has almost zero mean (e) and a very small standard de- 
viation (a). We obtain (e 512 ) = 10" 6 , (e 32 ) = 2 • 10~ 4 , 
(T512 = 2.3 ■ 10~ 4 , cr 32 = 3.71 ■ 10~ 3 . Although some pix- 
els around poles have more error, the fraction of these 



pixels on the sky for n s id = 32 is less than 0.5%. So 
for n s id c = 32, the error and standard deviation are suf- 
ficiently small, which enables one to calculate the a; m and 
to construct the cleaned map for KSP using Tegmark et al. 
method. 



3.2. Modified Tegmark et al. method 



We use Tegmark et al. method ( Tegmark et al.||2003 
et~aI1|200M [20081 |Gurzadyan et al.|J2009b 



to a 



eve 



Saha 



op a 

cleaned CMPTTCmap for eight WMAP bands: Ql, Q2, VI, 
V2, Wl, W2, W3, W4. This method, based on the power 
spectrum comparison, assigns a weight wj for each map i 
and multipole I. In order not to distort the original map, 



should obey the relation 



= i. 



(4) 



But a priori w\ could be any real number, including nega- 
tive. But KSP has to be within the interval < $ < 1 and 
the use of such negative weights to construct the cleaned 
ai m and then reconstruct the map could result in negative 
values of $ for some pixels. To avoid this problem, we nor- 
malize the weights by the following formula: 



(5) 



where w™ m is the minimal value of the original weights for 
fixed I multipole. It is easy to see that for the new weights 
w\ the relation J^. w\ — 1 is true. This is the unique linear 
relationship satisfying w\ > and minimally modifying the 
original w\. 

In Kolmogorov method, if a sequence of random num- 
bers T n obeys a theoretical distribution function F(x) then 
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Fig. 3. Cleaned K-map (top) and regions with higher KSP 
value ($ > 0.5) (bottom). The blue and green dots, outside 
Galactic a nd LMC regions, ar e poi nt sources li s ted in the 
catalogs of |Gold et al. ( 2011[ ) and Lanz et al. (20111, re- 
spectively. Seven of the blue dots are point sources list ed in 
the Planck Early Release Compact Source Catalogue ( Ade 
et al.||2011|). The remaining two red dots are still unidenti- 



fied possible point sources. 



N realizations of this sequence gives <&at, A at. The remark- 
able point of the method is that, <5>n has a uniform distri- 
bution and the mean value (<J > ) = 0.5. Therefore for a CMB 
K-map the value 0.5 is a natural threshold for distinguish- 
ing non-Gaussian areas from others, since Gaussian tem- 
perature pixels cannot have 4> > 0.5. Further, the cleaned 
K-map mean value and sigma are respectively ($) = 0.22, 
<7$ = 0.14, so pixels with $ > 0.5 mostly exceed the 3- 
cr region (see Fig. [5). In Fig. [3] one can see the Galactic 
disk, the Large Magellanic Cloud (LMC) and other pos- 
sible point sources. There are in total 398 non-Gaussian 
pixels among of 12288 (about 3.2%). Only 27 among them 
are outside the Galactic region (\b\ < 20) and the Large 
Magellanic Cloud (/ = 280.4136, b = 32.9310). So, these 27 
regions are point source ca ndidates. Indeed, most of them 

Gold et~aT1 ( |2011| ), |Lanz et al 



are fou nd i n the catalogs of 



(2011) and Ade et al. (201 



|. It is interesting that seven 



of the point sources are also listed in the Planck Early 
Release Compact Sourc e Catalogue as high reliability point 
sources (Ade et al. 2011 1. Two possible point sources remain 
unidentified in any known catalog. 

Another interesting feature of our cleaned K-map is 
that its distribution differs from the uncleaned maps in a 
minimal way. The restrictions from the weighting scheme 
Eq. Q keep the mean value almost unchanged. This means 
that structures in the cleaned map are preserved during the 
cleaning procedure, but for a lot of regions the KSP value 
is decreased. 




Fig. 4. Cosmic microwave background temperature and K- 
map correlation function power spectra. 

4. CMB and K-map power spectra 

Any full sky map can be repre sented via a series o f Legendre 
spherical functions Yi m (9, p) ( Hivon et al.||2002 ), where 

T(9,p) = ^ai, n Yi m (0,ip), 



l ,771 



air, 



T{9 7 ip) Y* m {6, ip) sin 9 d9 dp. 



(0) 



As easily shown, the coefficients Ci of the Legendre poly- 
nomials Pi(cos9) in the two point correlation function C(9) 
of the power spectrum are related to the a\ m by 



Ct = 



C{9) 



\ a lm a lm), 



l)CiPt(coa9). 



(7) 



l.rn 



To see if there is any correlation between the CMB tem- 
perature and the KSP parameter map, we degrade the res- 
olution of the CMB temperature map to n S id c = 32 and 
also normalize it to the same < T < 1 interval as 
So, instead of the original temperature we use a dimension- 
less parameter T, to prevent any discrepancy between tem- 
perature and Kolmogorov maps. For power spectr um and 
cross power spectrum calculations we use WMAP7 ( Jarosik 



et al. 



2011 ) eight different bands T and the Kolmogorov 
maps"(Ql, Q2, VI, V2, Wl, W2, W3 W4). Then, the fore- 
ground reduced map is obtained using Tegmark et al. ( 2003 ) 



weighting technique for triplets of different bands was used 
for the calculation of the power spectra both for the tem- 
perature and the KS P correlation function s, v ia standard 
techn i que as shown in Hinshaw et al. ( 2003 ) and Gurzadyan 
et al. ( 2009a I . The lines in Fig. [4] are the mean value of the 
Ci for different cross power spectra between different bands. 
The error bars are obtained as the variance square root of 
the different cross power spectra. 

The results for the foreground cleaned maps are shown 
in Fig. |4j It can be seen that the C; have different behavior 
and the main striking fact is that they are crossed at I ss 25. 

5. Cross power spectra of CMB temperature and 
K-maps 

For cross power spec tra calculation we us e the common 
technique described in Hinshaw et al. ( 2003 ) via taking the 
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frequency count for WMAP7 W KSP map 
frequency count for WMAP7 cleaned KSP map 
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Fig. 5. Histograms for WMAP7 W band and cleaned K- 
maps. 




90 100 



Fig. 6. Cross power spectrum of CMB WMAP7 T and $ 
maps. 



cross-correlation power spectrum coefficients for different 
type of ai m , as 



Ci — ( a lm a lm) 



(8) 



Here we again u se the foreground reduced maps (see 
Tegmark et al. 2003 ) for the calculation of the cross power 
spectra between T and since the original ones are very 
noisy, while the aims is to get the smallest possible error 
bars in the final power spectrum. The power spectrum es- 
timation is done without taking into account the Galactic 
disk plane region, i.e., we use the window function which 
is zero for the region 8 — ±20 for both T and $, and 
unity elsewhere. This cutting method influences mostly the 
even I inducing some unreasonable peaks (for example at 
I = 2, 12, 22, . . .) around t he window fun ction power spec- 
trum peaks. One could use Peebles (19731 method to reduce 
this effect on even I values and adjust odd I values but then 
one would have to calculate the power spectrum at least up 
to / = 250. The low resolution map of $ allows the accurate 
estimation of the power spectrum up to I — 97 which makes 
impossible the use of this method. Therefore we have to use 
only odd I values. 

In Figj6]no proper correlation can be seen between CMB 
WMAP7 T and $ maps. 



6. Conclusions 



Using Kolmogorov parameter for detection of point sources 
in CMB maps, we proceeded from the K-maps cleaned via 
the modified Tegmark et al. method. The novelty is that we 
applied that cleaning algorithm not to the usual tempera- 
ture maps, but for the K-maps. Other mask construction 
schemes are based on WMAP K and Kl band maps and are 
more affected by instrumental and other types of noises. It 
appears that for about 85% of the cleaned map Kolmogorov 
function has values in the interval 0.0 < $ < 0.2, which 
implies that the CMB maps are Gaussian with high pre- 
cision. The non-Gaussian pixels (with $ > 0.5) are rather 
rare (398), i.e., less than 4%. Of course this result is de- 
rived for the cleaned K-map but it implies that different 
types of noises add additional non-Gaussianities into the 
CMB maps, which may be anal yzed by the method s in[de| 
Oliveira-Costa et al.| ( |2004"| and |Rocha et al.| ( |2005[ ). These 
pixels when outside the Galactic disk region, |6| > 20, indi- 
cate the positions of point source candidates. While most 
of them have counterparts in existing catalogs, two of them 
are still unidentified. 



Another type of non-Gaussianity discussed here is the 
fluctuation of K-map. KSP is a statistical parameter so if 
one calculates KSP for fixed numbers of elements, then it 
must have statistical fluctuations. But for the same reason, 
the full sk y K-map should have un iform distribution which 
is not so (Gurzadyan et al. 20111. So in most cases these 
fluctuations have non-Gaussian nature. As it was shown in 



( |2009b[ ), the power spectrum for CMB 
K-map reveals some feature around Z = 25; the simplest 



Gurzadyan et al. 



possible explanation is the existence of a characteristic size 
of voids in the Universe, i.e., I — 25 corresponds to angular 
size 9 — 7.2°, as a result of the integrated Sachs-Wolfe 
effect. 



Some theoretical models dTakeuchi et al. 2012;|Komatsu 



1991[|Gangui|199 

from inflation era 



su 
4) 



|fc Spergel|2001||Salopek fc Bond 1990, 
predicts primordial non-Gaussianities 
This effect is very tiny but some authors tried to discover 
it through CMB bi-spectrum. We then used a new method 
to implement the cross-correlation between CMB temper- 
ature and K-map. Numerical modeling of such a problem 
was done in |Ghahramanyan et al. ( 2009 ) . It was shown that 
KSP is sensitive even to small departures from theoreti- 
cal distribution. Certain non-Gaussian perturbations would 
appear in K-map as KSP perturbations. Also, if one uses a 
proper theoretical distribution function, no correlation be- 
tween temperature and KSP should arise. Therefore in re- 
gions outside the Galactic disk certain correlation could ap- 
pear even in the presence of rather small non-Gaussianities. 
As one can see from Fig. [3] and Figj6j it is rather difficult to 
see any correlation, although the intersection of the correla- 
tion functions of T and $ a round / = 25 might be related t o 
certain symmetries, see e.g. Gurzadyan et al. (2008 2009a). 
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Table 1. Coordinates and ID-s of point sources detected 
in cleaned K-map. 



1 


b 


KSP value 


listed in 


source ID 


15.00 


58.92 


0.50250 


Planck 


- 


43.59 


27.28 


0.50681 


- 


- 


45.00 


49.70 


0.50025 


WMAP7 




63.28 


41.81 


0.57482 


WMAP7 


GB6 J1625+4134 


67.50 


-34.23 


0.52509 


Lanz et al. 


GB6 J2148+0657 


72.69 


70.91 


0.54416 


WMAP7 


GB6 J1419+3822 


85.78 


-38.68 


0.83258 


WMAP7 


GB6 J2253+1608 


97.03 


24.62 


0.56636 


WMAP7 


GB6 J1841+6718 


99.84 


38.68 


0.51425 


WMAP7 


J1659+6827 G 


99.84 


24.62 


0.56640 


WMAP7 


J1842+6808 QSO 


105.47 


24.62 


0.52249 


WMAP7 


GB6 J1927+7357 


106.50 


44.99 


0.61560 


Planck 


- b 


132.19 


66.44 


0.50769 


Planck 


- h 


157.50 


-20.74 


0.54468 


Lanz et al. 


GB6 J0336+3218 


158.91 


-22.02 


0.54752 


Planck 


- 


195.47 


-32.80 


0.71974 


WMAP7 


PMN J0423-0120 


210.94 


-20.74 


0.70129 


Planck 


- b 


213.75 


-20.74 


0.62933 


Planck 


- b 


238.33 


-49.70 


0.54351 


WMAP7 


PMN J0406-3826 


246.09 


24.62 


0.51170 


Planck 


b 


251.72 


-32.80 


0.59206 


WMAP7 


PMN J0526-4830 


258.75 


-72.39 


0.65437 






277.03 


-35.69 


0.50510 


WMAP7 


PMN J0537-6620 


282.27 


73.87 


0.51290 


WMAP7 


GB6 J1230+1223 c 


283.50 


75.34 


0.73977 


WMAP7 


GB6 J1230+1223 c 


288.53 


64.95 


0.83234 


WMAP7 


QSO J1229+0203 


304.77 


57.40 


0.78074 


WMAP7 


PMN J1256-0547 


a unidentified point sources 







possibl e point sou r ce ide ntified in Planck Early Release Sources 
catalog |Ade et al.| < |2011[ 
c unidentified in WMAP7 catalog 
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